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IvsCeLaSit?Cchr^e  SCheTOeS  for  vasFiou  1-D  hyperbolized  dissipative 
systems  with  stiff  source  terms  w?re  developed  and  tested  with  success, 

Pr€COn^tionill9  matrix  that  maintains  the  largest  possible 

enti  Mach-11  KS  eigenvectors  of  the  preconditioned  system  for  the 
enti_e  Mach  number  range,  was  developed  in  order  to  prevent 

the  observed  stagnation-point  instability,  and  tested  with  success. 

A  Navier-Stokes  preconditioning  matrix  that  remains  stable  and 

a^fr  rV°r  JU  MaCh  nurobers  and  Reynolds  numbers  was  developed 
and  tested  with  success. 
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Final  Technical  Report 

Explicit  integration  schemes  for  the  hyperbolized 
Navier-Stokes  equations 

Grant  No.  F49620>93-1-0417 

Principal  Investigator: 

Bram  van  Leer 

Department  of  Aerospace  Engineering, 

University  of  Michigan 

27  March  1997 


1  Objectives 


Tins  project  considered  the  application  of  local  preconditioning  to  a  stiff  hyperbolic 
system  describing  viscous  conducting  flow  with  a  finite  relaxation  time.  It  was  funded 
y  an  Augmentation  Grant  tacked  on  to  the  larger  project  “Local  preconditioning  of 
he  Euler  equations  and  its  numerical  applications,”  funded  under  Air  Force  Grant 
No.  F49620-92- J-0 1 58- DEF .  During  the  second  funding  year,  part  of  the  effort  in  the 
augmentation  project  was  diverted  to  the  larger  subject  of  preconditioning  for  the 
Euler  and  Navier-Stokes  equations,  in  order  to  help  solve  the  problem  of  stagnation- 
point  instability,  which  was  halting  any  progress  in  the  application  of  local  precondi- 
tionmg  This  situation  persisted  during  the  third  year  and  following  no-cost  extension 
peno  (  / 1-12/31/97),  particularly  because  another  preconditioning  instability  one 
lor  very  low  cell  Reynolds-numbers,  was  encountered.  The  reformulated  goals  for  the 
augmentation  project  therefore  were: 


1.  To  develop  accurate,  fast  and  robust  integrators  for  hyperbolized  dissipative 
systems,  but  without  the  use  of  local  preconditioning; 

2.  To  develop  a  local  preconditioning  for  the  Euler  equations  with  non-degenerate 
eigenvector  system  for  the  entire  Mach-numbcr  range  (the  degeneracy  causes 
the  stagnation-point  instability); 

3.  1©  develop  a  local  preconditioning  for  the  original  (non-hyperbobzed)  Navier- 
.tokes  equations  that  is  effective  and  stable  for  all  combinations  of  Mach  and 
ceil  Reynolds-numbers. 
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2  Final  status  vs.  goals 


oa)  1:  The  development  and  testing  of  integration  schemes  for  hyperbolized  dis- 
Vf*  1  Ve,  systems  with,  stiff  source  terms  was  rounded  off  with  the  production  of  the 
CTa  theS1S  °fMoh,t  Arora>  defended  in  February  1996.  The  systems  studied  were: 

(a)  1-D  hyperbolized  heat  equation  (2  eqs.); 

(b)  hyperbolized  Burgers  equation  (2  eqs.); 

(c)  1-D  Broad  well  gas  (3  eqs.); 

(d)  1-D  Euler  equations  with  heat  conduction  (3  eqs.); 

(e)  1-D  10-momcnt  equations,  equivalent  to  N-S  without  heat  conduction  (1  eqs.). 

W„?!  2:  Afn?,1‘riX  1'V!l0p'!d  preV,0US‘y  (md'r  *k«  parent  grant)  by  Dohyung  Lee 
.  ucces.sfii.lly  applied  to  a  standard  stagnation-point  test  flow.  Convergence  ac- 

celeration  was  obtained  where  the  Van  Leer-Lec-Roe  and  Turtel  matrices  produce  a 

mofal  TqT  ‘frfoor?’  TUs  “  Very  Promisin«-  si““  ‘he  “fix”  proposed  by  Dar- 
”  f,  '  °d  Schmid  (1995)  requires  the  introduction  of  a  rather  arbitrary  lower  limit 

r  e  Mach  number  which  significantly  complicates  the  algebra  in  designing  the 
couesponding  artlhcial-viscosity  matrix.  The  Lee  matrix  has  a  non-degenerate  (non 

LTi  mcdTfr0rtrCtUremtr  "  X  °’  b“‘  “°l  f“'  M  T  li  «he  Van  Leer  matrix 
funct  on  f  y  2®  t W  ?y-  TbiS  SUSgCStS  »  blrad  two  matrices  as  a 

fo  1  M  6  WaCh  MmW‘  A  matrix  was  foulld  ‘hut  »  close  to  Lee's  matrix 

¥  close' to  T  P°!!  “ve  drfmf  -Le*  s  is  not)  and  approaches  Van  Leer’s  matrix  for 
close  to  1  with  improved  eigenvector  structure  for  the  entire  range  of  M  This 

mac  to t T  -- Itchtained  in  the  project,  and  has  not  yeAeen  tested  in 

rNmernbert. 96  t,0n  «“  be  f°’md  D'  L"’s  d^nded 

®0al.  3:.  F0l!°"'mS  1  S'JS?“t,on  b-v  S-  Veniateswaran  (NASA  ARC.  1995.  private 

to rCS1  "  T  T  "T  01 1  Simplificd  Ni“‘Cr  St0k«  analysis 

to  put  Reynolds-number  dependence  into  the  Van  Leer  preconditioner.  This  does 

equahac  the  various  time  scales  in  the  N-S  eqnations,  but,  unfortunately,  produces  an 

unstable  mode  for  very  low  cell  Reynolds-numbers.  D.  Lee  found  that  the  precondi- 

tfon  formTtl  C1  -V  addinS  *°  “  tbe  matdX  “efficients  of  the  leading  dissipa. 

on  terms  that  would  also  arise  in  point-Jacobi  relaxation.  This  preconditioning  was 

successfuHy  .mplemcnted  in  the  computation  of  a  flat-plate  laminar  boundary-laver 

laver  t! ^  “eU-^Pect-rntio  going  up  to  1700  at  the  plate,  and  a  turbulent  boundary 
layer.  These  results  are  contained  in  D.  Lee’s  thesis.  ^ 
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3  Project  accomplishments/highlights 

Ths  research  effort  spent  in  the  3.5  years  of  this  augmentation  project  has  produced 
a  significant  improvement  of  local  preconditioning  tor  the  standard  NavZ-Stokes 
e,ua  ions  A  matrix  was  developed  that  equalises  the  tune-scales  imp.  "  in  the 

with  -trallr  “  ‘TilCh  rd  Cell‘R'y,loldii  numb«s.  and  leads  to  stable  calculations 

hln  f .  h^i  aCMl'ra'ed  “^gence.  Such  a  matrix  produces  not  only  a  single-grid 

current.  N  _al%an  mdep(!nd“*  multi-grid  benefit.  This  is  highly  desirable,  because 

stiff  a  do,:TS  ^  C°dCS  are  °°t0n0US  for  th"r  slow  convergence.  Robustness  is 
a  point  of  concern;  more  testing  is  needed* 

“belh-ht-  “T  ?  Eul'r  pr'c01“iiti0n“S  we  may  have  found  a  matrix  with  a  sufficiently 
healthy  effect  on  eigenvector  structure  so  as  to  yield  stable  calculations  even  in  the 

be  mldV  d‘  sli'gMtlon  Pomt-  More  ‘“ting  is  still  needed;  the  matrix  will  eventually 

calculations.  “  Reynolds-“““ber  dependence,  for  use  in  Navie, -Stokes 

Z  T°  f  hypcrboli;!'d  dissipativ*  systems,  the  integration  method  developed 
was  Znd  WhW°r  VT  ,  *  "ide  VariCty  °f  1-D  Systcros  'Stations  and 

1  i  Z  '  ,S1‘k'TaClory  K  ‘“'•curate  and  fast,  and  extremely  robust  In  the 
test  calcinations  it  handled  relaxation-time-scales  ranging  over  a  factor  10“  The 

mam  challenge  is  to  extend  the  approach  to  multidimensional  systems.  The  benefit 
of  the  hyperbolized  approach  is  that  infinite  signal  speeds  (as  in  the  standard  diffu- 
to  eo'ThT  4re  a™lded.  relieving  the  severe  time-step  restriction,  and  the  need 
theSsou  °b1  y  lmpl,Clt- . A  ,ocally  implidt  integration  technique  then  takes  care  of 
r*0”'  ThiS  may  be  the  Way  °f  the  **«•  of  computing  any 

“d  «ow  at 


4  Personnel  supported 

Brain  van  Leer  (PI),  professor, 

Philip  L.  Roe,  professor, 

John  F.  Lynmn  doctoral  candidate  (defense  17  May  1995) 
Dohyung  Lee,  doctoral  candidate  (defense  November  20,  1996 
ohit  Arora,  doctoral  candidate  (defense  February  8,  1996) 
4.11  in  the  same  department. 


5  Publications 

1  Al°til  fd  PhiIip  L  R°e-  ” Characteristic-based  numerical  algorithms 

for  stiff  hyperbolic  relaxation  systems,"  presented  at  the  15th  International 
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oZtToZ  °P  Num,ericaI  ^ethods  in  Fluid  Dynamics,  Monterey,  CA,  June 
24-28,  1996.  Proceedings  to  be  published  with  Springer. 

2.  Mohit  Arora,  “Explicit  characteristic-based  high-resolution  algorithms  for  hv- 
MK^an!O1^7atl0n  kWS  With  Stiffs0UrCe  termS’”  Ph-D'  ^  University  of 

3  M°7)  ^Strategy  for  -hyperbolic  systems  with  stiff  source  terms,”  pre- 

sented  (by  R  L.  Roe)  at  the  ICASE/LaRC  workshop  on  Barriers  and  Challenges 

Kluwer  mpl°n’  VA)  AugUSt  5*7’  1996'  Proceedings  to  be  published  with 

4’  anJ'rLyTpaild  BJ-Ian  Leei’  “MuItjgfid  Euler  Solutions  with  Semi-Coarsening 
and  Loca  Preconditioning’ ,  in:  Proceedings  of  the  14th  International  Confer¬ 
ence  on  Numerical  Fluid  Dynamics,  Springer,  1995. 

5.  J.  F.  Lynn  and  B.  van  Leer,  “A  Semi-Coarsened  Multigrid  Solver  for  the  Euler 

qUati°nS  Witt  Local  P^conditioning”,  AIAA  Paper  95- 

6.  John  F.  Lynn,  Bram  van  Leer  and  Dohyung  Lee,  “Multigrid  solution  of  the 
r  ;  equations  with  local  preconditioning,”  presented  at  the  15th  International 

S0niOQrnCp°n  Nl;mericaVMethods  m  F,uid  Dynamics,  Monterey,  CA,  June  24- 
-8,  199G.  Proceedings  to  be  published  with  Springer  Verlag. 

D.  Lee  “Local  Preconditioning  of  the  Euler  and  Navier-StokesEquations Ph  D 
thesis,  university  of  Michigan,  1996. 


6  Interactions  /  transitions 

A.  Participation/presentations,  etc. 

The  Principal  Investigator  Bram  van  Leer  gave  numerous  CFD  seminars  and 
other  presentations  (without  ensuing  publications)  in  which  results  on  precon- 

•  W6v  prese^led  u  Starting  with  the  fall  of  1995,  such  presentations 

rnduded  positive  results  obtained  in  the  augmentation  project.  These  are  listed 

•  September  16,  1995,  “Local  preconditioning  of  the  Euler  and  Navier-Stokes 
equations,  invited  lecture  at  the  Conference  on  Numerical  Methods  for 
the  Euler  and  Navier-Stokes  Equations,  Montreal,  Quebec,  Canada 

•  1996;USAF  Philips  Lab.,  Kirtland  AFB,  Albuquerque,  NM, 
CrD  in  the  nineties  and  beyond’5 

•  March  21,  1996,  CFD  Research  Corp.,  Huntsville,  AL,  “Is  CFD  dead?” 
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•  March  28,  1996,  Texas  A&M,  College  Station,  TX,  “Is  CFD  dead?” 

•  April  24,  1996,  MacNeal-Schwendler  BV,  Gouda,  Netherlands,  “Is  CFD 
dead?’ 

•  April  29,  1996,  “Does  industry  need  fanciful  CFD  methods?**  plenary  lec- 
tiu^  at  the  BRITE/EURAM  Workshop  on  Multi-Dimensional  and  Sparse- 
Gnd  Methods,  Von  Karman  Institute,  St.  Genesius-Rode,  Belgium. 

•  May  2,  1996,  MacNeal-Schwendler  BV,  Gouda,  Netherlands,  “Local  pre¬ 
conditioning  tor  the  Euler  and  Navier- Stokes  equatons.” 

•  May  6,  1996,  CWI  (Center  for  Mathematics  and  Informatics),  Amsterdam, 
Netherlands,  “Numerical  integration  of  stiff  hyperbolized  systems:  whv 
and  how.” 

•  August  7,  1996,  lowards  the  ultimate  local  preconditioning  matrix:  Is  the 
end  in  sight?”  Invited  lecture  at  the  ICASE/LaRC  workshop  on  Barriers 
and  Challenges  in  CFD,  Hampton,  VA,  August  5-7,  1996. 

B.  Consultative  and  advisory  functions,  etc. 

•  Bram  van  Leer  visited  the  USAF  Phillips  Lab.,  Kirtland  AFB,  Albu¬ 
querque,  NM,  on  February  15.  Discussions  included  the  use  and  numerical 
treatment  of  hyperbolized  dissipative  systems, 

•  Bram  van  Leer  is  a  regular  consultant  to  ICASE  (Institute  for  Computer 

Applications  m  Science  and  Industry),  located  at  NASA  Langley  Research 
Center,  and  has  advised  NASA  branches  extensively  on  the  uses  and  im¬ 
plementation  of  local  preconditioning  of  the  Euler  and  Navier-Stokes  equa¬ 
tions.  n 

C.  Transitions 

Di.  Veer  Vatsa  (NASA  Langley  Research  Center)  has  successfully  implemented 
m  the  aerodynamics-analysis  code  TLNS3D  a  preconditioner  of  the  family  that 
includes  the  matrices  of  D.  Lee  and  B.  van  Leer,  which  remain  stable  in  the  pres¬ 
ence  of  a  stagnation  point.  Implementation  of  these  new  matrices  is  scheduled 
to  take  place  during  a  collaborative  project  funded  under  Advanced  Subsonic 
technology,  during  the  summer  of  1997. 


7  New  discoveries,  etc. 

None  patentable. 


8  Honors  and  awards 

•  Bram  van  Leer  was  made  Fellow  of  the  A1AA  in  May  1995. 
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•  Brain  van  Leer  received  the  UMI  College  of  Engineering’s  Award  for  Excellence 
in  Research,  January  1996. 

•  Phil  Roe  was  made  Fellow  of  the  AIAA  in  May  1996. 
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